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I. INTRODUCTION 


A comparison of QCD theoretical results with experimental data is often based on the 
concept of quark-hadron duality, which establishes a bridge between quarks and gluons, a 
language of theoreticians, and real measurements with hadrons performed by experimental¬ 
ists. The idea of quark-hadron duality was formulated in the paper by Poggio, Quinn, and 
Weinberg Q] as follows: Inclusive hadronic cross sections, once they are appropriately aver¬ 
aged over an energy interval, must approximately coincide with the corresponding quantities 
derived from the quark-gluon picture. For many physical quantities and functions the corre¬ 
sponding interval of integration involves an infrared region and in this case nonperturbative 
effects may play an important role in their description. 

In this paper we consider the following quantities and functions. 


The ratio of hadronic to leptonic r-decay widths in the vector channel; 


ds 


M3 


fit = / tTj I-TF? 1 + T7J 


M3 


2s 


M3 


( 1 ) 


the so-called “light” Adler function, which is constructed from r-decay data 

m . 


Q: 


D{Q^) = -Q" = Q^ I ds 


dQ^ 


(5 + Q2)2’ 


( 2 ) 


the smeared function 


ion [l|: 


OO 

A C 

-Ra(’S) = — ds' 


R{s') 


(s — ’ 


( 3 ) 


the hadronic contribution to the anomalous magnetic moment of the muon (in the 
leading order in electromagnetic coupling constant): 


OO 


( 4 ) 


where K{s) is the vacuum polarization factor given by (ITH|) below; 


• and the strong interaction contribution to the running of the hne structure constant: 




a(0) 

Svr 


M|P 



Rjs) 
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A common feature of all these quantities and functions is that they are dehned through the 
function R{s), the normalized hadronic cross-section, integrated with some other functions. 
By dehnition, all these quantities and functions include an infrared region as a part of the 
interval of integration and, therefore, they cannot be directly calculated within perturbative 
quantum chromodynamics (pQCD). 


The method that we use here to describe the quantities and 
is based on the analytic approach to QCD suggested in 
allows one to describe self-consistently the timelike region 
integration in Eqs. O-dSD- It incorporates the required analytic properties and leads to an 
integral representation for R{s). We formulate a model that also incorporates a summation 



unctions mentioned above 
The analytic approach 
, which is represented in the 


of threshold singularities 
light quark masses. 


Q 


and takes into account the nonperturbative character of the 


II. METHOD AND BASIC RELATIONS 


A. Analytic perturbation theory 

Analytic perturbation theory (APT) jsl is based on the analytic approach to QCD. In this 
approach, in contrast to the behavior of the perturbative running coupling, the Euclidean 
analytic coupling has no unphysical singularities. The ghost pole and corresponding branch 
points, which appear in higher orders, are absent. APT preserves the correct analytic 


properties of such important objects as the two-point correlation 


a well-dehned algorithm for calculating higher-loop corrections joj. In APT, processes with 
typical spacelike and timelike momenta are described self-consistently lO; lld| and, for 
example, inclusive r-decay can be described equivalently either in terms of Minkowskian or 
Euclidean variables j^. In the framework of APT, the theoretical ambiguity associated with 
the choice of renormalization scheme is dramatically reduced. 

In the APT scheme the QCD contributions d{z) and r(s) to the functions D oc 1 + d 
and R (X 1 + r, respectively, are expressed in terms of the effective spectral function p(cr) as 


unction and also provides 


d{z) = - 


da 


TT 


a 


p{a ), r(s) = - [ —p{(^) ■ 

TT a 


( 6 ) 


The APT spectral function p{a) is defined as the imaginary part of the perturbative function 
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dpt{z) and in the third order can be written in the form p{a) = Po(o') + di^i(cr) + 
where Pk(o') = Im[ap^^(cr + ie)], a = cTs/tt, and dk are the coefficients of the perturbative 
expansion of the D-function. 

The function po(<^) dehnes the analytic spacelike, A{z), and timelike, 2l(s), running 
couplings as follows 


A(z) = 


TT 


da 


a — z 


go{a ), 2t(s) = 


TT 


da 


a 


Qo{a). 


(7) 


As has been argued from general principles in [l]|, the behavior of these couplings cannot be 
the same, i.e., they cannot be symmetrical in the spacelike and timelike domains, A{—z) ^ 


2l(z). The analytic and perturbative couplings have been compared in |^. In analyzing 
hadronic processes with characteristic spacelike and timelike momenta, it is necessary to 
take into account this lack of symmetry between the behavior of the running coupling in the 
Euclidean and Minkowskian regions. 

The analytic running coupling has no unphysical singularities and possesses the correct 
analytic properties, arising from Kallen-Lehm^n^nalyticity reflecting the general principles 
of the theory The one-loop APT result is 


A^^\z) = ApAz) + 


4 A2 
PqA'^ + z' 




Po 


1 1 ln(s/A2) 

-arctan- 

2 TT TT 


( 8 ) 


where apt^(z) = ds{z)/'n' = 4/ [Po\n{—z/AP)] and /3o = 11 — 2//3 is the first coefficient of 
the renormalization group /d-function. 

Both the couplings m have the same infrared fixed point ^(0) = 21(0) = 4//3o. This value 
is dehned by the leading contribution (jH)) and is not altered by higher-order corrections. The 
regular behavior in the infrared region of (z) is provided by the power term in (jH)) which 
is invisible in the perturbative expansion. The reason for the regularity of the coupling 
2lW(s) is connected with the summation of 
role in analyzing various hadronic processes 


the so-called tt -terms that play an important 


iaiiaiiaiij,ii5i,ii6|. 


B. Resummation of threshold singularities 

In describing a charged particle-antiparticle system near threshold, it is well known from 
QED that the so-called Coulomb resummation factor plays an important role. This resum¬ 
mation, performed on the basis of the nonrelativistic Schrodinger equation with the Coulomb 
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potential V{r) = —a/r^ leads to the Sommerfeld-Sakharov S'-factor HQ . In the thresh¬ 
old region one cannot truncate the perturbative series and the S'-factor should be taken into 
account in its entirety. The S'-factor appears in the parametrization of the imaginary part 
of the quark current correlator, which can be approximated by the Bethe-Salpeter ampli¬ 
tude of the two charged particles, Xbs( 3^ = 0) 19|. The nonrelativistic replacement of this 
amplitude by the wave function, which obeys the Schrodinger equation with the Coulomb 
potential, leads to the appearance of the resummation factor in the parametrization of the 
i?(s)-function discussed above. 

For a systematic relativistic analysis of quark-antiquark systems, it is essential from the 
very beginning to have a relativistic generalization of the S'-factor. A new form for this 
relativistic factor in the case of QCD has been proposed in 


^(x) = 


^(x) 


^(x) = 


TT a 


(9) 


1 - exp [-X(x)] ’ sinhy’ 

where x is the rapidity which related to s by 2m cosh x = \/s, « —> 4 q!s/ 3 in QCD. The 
function X{x) can be expressed in terms of u = — Am?/s\ X{x) = 7ra-\/l — v‘^/v. 

The relativistic resummation factor m reproduces both the expected nonrelativistic and 
ultrarelativistic limits and corresponds to a QCD-like Coulomb potential. Here we consider 
the vector channel for which a threshold resummation S'-factor for the s-wave states is used. 
For the axial-vector channel the P-factor is required. The corresponding relativistic factor 
has recently been found in Q 

To incorporate the quark mass effects one usually uses the approximate expression pro¬ 


posed 


in i 


21| above the quark-antiquark threshold 


7J(s) = T(v) [1 + 9 (i')i'(s)] , 


( 10 ) 


where 


T{v) = V- 


4:71 

9(1.) = ^ 


71 


3 + V /' TV 
2 


3 

4tt 


Vf = \^ 


4m?j 


( 11 ) 


The function g{v) is taken in the Schwinger approximation 221. 

One cannot directly use the perturbative expression for r(s) in Eq. (uni), which contains 
unphysical singularities, to calculate, for example, the Adler P-function. Instead, one can 
use the APT representation for r(s). The explicit three-loop form for rAPT(s) can be found 
in j^. Besides this replacement, one has to modify the expression m in such a way as to 
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take into account summation of an arbitrary number of threshold singularities. Including 


the t 

(HUD 


ireshold resummation factor Q leads to the following modihcation of the expression 
[ 1^ for a particular quark flavour / 


TZfis) = [Roj{s) + Rij{s)] 0(s - 4m}), 


( 12 ) 


Ro{s)=T{v)S{x), Ri{s)=T{v) 


rAPT{s)g{v) - 2^(x) 


The usage of the resummation factor © reflects the assumption that the coupling is taken 
in the V renormalization scheme. To avoid double counting, the function Ri contains the 
subtraction of X(y;). The potential term corresponding to the Rq function gives the principal 
contribution to TZ(s), the correction Ri amounting to less than twenty percent for the whole 
energy interval 


C. Quark masses 


The following considerations suggest the behavior of the mass function of the light quarks 
in the infrared region. A solution of the Schwinger-Dyson equations 24, 3,1^1 demonstrates 


a hxed infrared behavior of the invariant charge and the quark mass function. The mass 
function of the light quarks at small momentum looks like a plateau with a height approx¬ 
imately equal to the constituent mass, then with increasing momentum the mass function 
rapidly decreases and approaches the small current mass. 

This behavior can be understood by using the concept of the dynamical quark mass. This 
mass has an essentially nonperturbative nature. Its connection with the quark condensate 
has been established in j^|. By using ^analysis based on the Schwinger-Dyson equations 


a similar relation has been found in 


It has been demonstrated in 


mass-shall one has a gauge-independent result for the dynamical mass 


= —7ra5(0|gg|0). 


that on the 


(13) 


A result obtained in Q demonstrates the step-like behaviour of the mass function. The 
height m of the plateau is given by the quark condensate (I13|) . According to these results 
it is reasonable to assume that at small the function m{p^) is rather smooth (nearly 
constant). In the region > 1-2 GeV the principal behavior of the function m{p^) is 
dehned by perturbation theory with the renormalization group improvement. 
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FIG. 1; Effective quark mass. 


TABLE I: Typical values of and M^. 


/ 

u 

d 

s 

c 

b 

t 

ml (GeV) 

0.004 

0.007 

0.130 

1.35 

4.4 

174 

Mi (GeV) 

0.260 

0.260 

0.450 

1.35 

4.4 

174.0 


The following analysis was performed by nsing the model mass fnnction m{p^) that is 
shown in Fig. ^ We take the cnrve that connects the points Pa and Pb to have the form 
^3/ ^p 2 _ j^ 2 y parameters mo are taken from the known valnes of the rnnning (cnrrent) 
masses at pf, = 2 GeV. The qnantities considered here are not too sensitive to the parameters 
of the heavy qnarks and we take for c, b and t qnarks [p^) = mg = Mg = const. The 
valnes of mg at 2 GeV and typical valnes of Mg are shown in Table Q] 

III. PHYSICAL QUANTITIES AND FUNCTIONS GENERATED BY R{s) 

In this section we apply the model we have formnlated to describe the physical qnantities 
and fnnctions connected with i?(s), described in the Introdnction. 

A. Inclusive decay of the r-lepton 

The ratio of hadronic to leptonic r-decay widths in the vector channel is expressed by 
Eq. (ID), where = 3 |I4dP *S'ew/2, |I4d| = 0.9752 ± 0.0007 is the GKM matrix element, 
Sew = 1.0194 ± 0.0040 is the electroweak factor, and Mr = 1776.99lo;26 MeV is the mass of 
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the r-lepton 31 
for this ratio is 


.• T 

3 , 


re experimental data obtained by the ALEPH and OPAL collaborations 


33, 


^ALEPH _ ^ 775 ^ Q_Q77^ t^opal _ 7 7g4 ^ q 


In our analysis we use the nonstrange vector channel spectral function obtained by the 
ALEPH collaboration j3l and keep in all further calculations the value as the 

normalization point. The range of estimates are obtained by varying the quark masses in 
the interval Mg= 260 ± 10 MeV (this band is hxed rather dehnitely by the H-function 
considered below) and Mg = 450 ± 100 MeV. The results for are given below. 


B. Hyfunction 

The experimental information obtained by the ALEPH and OPAL collaborations allows 
us to construct the nonstrange vector channel “experimental” H-function. Within the an¬ 
alytic approach this function has been analysed in P|. Here we improve our method of 
constructing the “light” H-function by taking into account the global duality relation. We 
demonstrate that this Euclidean object is useful from the point of view of dehning the 
effective masses of the light quarks. 

In order to construct the Euclidean H-function o we use for R{s) the following expression 


R{s) = e{so - s) + e{s - So) • 

The continuum threshold sq we hnd from the global duality relation 


(14) 


So So 

0 0 


theor/ 


s). 


This gives sq — 1.5 GeV^. The value of sq agrees with the results of papers 


similar value of the continuum parameter is used in the QCD sum rules 


39 


(15) 

^ 0,0- A 

^|. Note, for 


some parameters there are two possible solutions of the duality condition m- We exclude 
the second solution, sq — 2.5 GeV^, at this stage of the analysis due to the requirement of 
describing, in a self-consistent manner, different experimental data. 

The low energy r-data in the nonstrange vector channel results in the curve for D{Q^) in 
Fig. El In this hgure we also plot three theoretical curves corresponding to masses of the light 
quarks of 150, 260 and 350 MeV. Fig. El demonstrates that the shape of the infrared tail of the 
H-function is quite sensitive to the value of the light quark masses. Note the experimental 
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FIG. 2: D-function for m = const. 


FIG. 3: L)-function for m = m{p^). 


D-function turns out to be a smooth function without any trace of resonance structure. 
The D-function obtained in Ref. from the data for electron-positron annihilation into 
hadrons also has a similar property. 


1.2 


0.8 


0.6 


0.4 


0.2 



0.35 


FIG. 4; F)-function surface at Qq = 0.5 GeV vs. parameters Mq and pa- The plane corresponds 
to h)expt(Qo)' 


The measured quantity dehned in Eq. m is less sensitive to and rrid values than 
the infrared tail of the Dy-function. Varying the light quark masses over a wide range one 
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finds = 1.79 for = 150 MeV and = 1.66 for 350 MeV. The values of 


masses 


iriu = TTid ~ 260 MeV agree with the experimental value R^ = 1.775±0.017 m . The values 
of the light quark masses are close to the constituent quark masses and therefore incorporate 


nonperturbative effects. T 
the analysis performed in 


hese va. 


36 


45 


ues are consistent with other results 


42 


43 


44l | and with 


46 \. 


A result for the H-function that is obtained by using the mass function m{p^) with 
parameters dehned in Tableland pa = 0.8 GeV is shown in Fig. El Thus we obtain results 
that are rather close to the results obtained for m(p^) = const = 260 MeV. In Fig. |3]we 
plot a 3-dimensional graph of the function D{Q‘^) as function of the parameters Mq and 
a. The plane corresponds to the experimental value i7expt(Qo) ~ at Qq = 0.5 GeV. 
Fig. m demonstrates that for pa > 0.4-0.5 GeV the curve of intersection of the surface 
Mo, Pa) and the plane is approximately a straight line, corresponding to Mq = 260 
MeV. The large Pa-limit reproduces the results with m{p^) = const. 


C. Smeared i^A-function 

To compare experimental and theoretical results from the point of view of the quark- 
hadron duality, in [l| it was proposed to use the smeared function -Ra(s). Instead of the 
Drell ratio R{s) dehned in terms of the discontinuity of the correlation function n(g^) across 
the physical cut 

7?(s) = ^ [n(s ± ie) - n(s - ie)], (16) 

the smeared function i?A(s) is dehned as 

Rii{s) = ^ [n(s±*A)-n(s-zA)], (17) 

with a hnite value of A to keep away from the cut. If A is sufficiently large and both the 
experimental data and the theory prediction are smeared, it is possible to compare theory 
with experiment. 

Equation (HZD and the dispersion relation for the correlator n(g^) give the representation 
m- Note that the smeared function Ra{s) is dehned both in the Minkowskian region of 
positive s, where a trace of resonances still remains for not too large A, and in the Euclidean 
domain of negative argument s, where like the Adler function D{Q‘^) the function R/\{s) is 
smooth and monotone. 
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FIG. 5: Smeared function for A = 0.5 GeV^. 



FIG. 6: Smeared function for A = 1.0 GeV^. 


As with the Adler function we will construct the “light” experimental function R^{s). 
For this purpose we match the experimental data taken with s < sq to the theoretical result 
taken with s > sq as in m- The value Sq — 1.6 GeV^ is found from the duality relation 

m- 

For the charm region the value of A is about 3 GeV^. An adequate choice in the case 
of the light smeared function is A ~ 0.5-1.0 GeV^. In Figs. El and IHl the experimental and 
theoretical curves for A = 0.5 GeV^, A = 1.0 GeV^ and m = m{p^) are shown. Let us 
emphasize that, for reasonable values of A, in the spacelike region (s < 0) there is a good 
agreement between data and theory starting from s = 0. 


D. Hadronic contribution to a,, 


The hadronic contribution to the anomalous magnetic moment of the muon in the leading 


order in the electromagnetic coupling constan 

137.035 99911(46) j^|, and (see, for example, 

1 

K{s) = [ dx 


; is dehned by (jH), where a ^ = Q!(0) ^ = 


m 


a;^(l — x) 


+ (1 — x)slm‘i 

0 

The muon mass is = 105.7 MeV. 

The expression (0 can be rewritten in terms of the H-function 

1 


(18) 


„had _ 1 1 


dx 


(1 — x){2 — x) D 


X 9 

-m;, 

1- X ^ 


(19) 


It is should be emphasized that the expressions 0 and (Hi are equivalent due to the 
analytic properties of the function n(g^). If one uses a method that does not maintain the 
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required properties of expressions (jj)) and (fTT?|l will no longer be equivalent and will 

imply different results (see for details). This situation is similar to that which occurs 
in the analysis of inclusive r-decay where the initial integral, performed over an interval 
including a nonperturbative region, for which a perturbative QCD calculation is not valid, is 
transformed based on the analytic properties into a contour representation. Within APT one 
is justihed in doing this, and can use equally well either the expression (0 or the expression 

P|). 


TABLE II: Dependence of on light quark masses. 



Wd ^ ;^Ql0 

niq (MeV) 

LO 

NNLO 

q = u,d q = s 

mg = const mg 

/ const 

mg = const mg / const 

250 400 

736 

760 

725 763 

250 500 

716 

736 

705 726 

260 400 

691 

715 

682 711 

260 500 

671 

690 

661 685 


The value of is not very sensitive to the values of the heavy quark masses, which 
we take as given in Table IH The relative contributions of u and d quarks are about 72 and 
19 %, respectively. The relative factor of 4 between u and d contributions is explained by 
the ratio of quark charges. The relative contribution of the s-quark to is about 5-9 % 
for Mq = 400-500 MeV. The contribution of the c-quark is about 2 %. Contributions of b 
and t quarks are very small. 

There is a signihcant dependence on the mass parameters of the light quarks. This de¬ 
pendence we illustrate in Table E] In our calculations we take into account the matching 
conditions at quark thresholds according to the procedure described in j^. The mass pa¬ 
rameters of u and d quarks are hxed rather well by the infrared tail of the light D-function 
and the value of . If we take for the parameter Mq’’^ in the function m = m{p^) the best 
£t value 260 MeV and vary Mg = 400-500 MeV, we get 


,had 


= (698 ± 13) X 10"^^ (20) 

Alternative “theoretical” values of are extracted from annihilation and r decay 
data: (696.3 ± 6.2exp ± 3.6i.ad) x 10“^° (e'''e“-based) 2, which is 1.9a below the BNL 
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H; 


experiment j^J; (711.0±5.0exp±0.8rad±2.8s'f/(2)) x 10 ('r-bas ed') J48j |. which is within 0.7a 
of experiment; and (693.4±5.3exp±3.5rad) x 10“^° (e'''e“-based) j^, 2.7a below experiment. 


An even lower value (692.4 ± 5.9exp ± 2.4i.ad) x 10“^° is given by [^. The quantity 
is rather sensitive to the light quark mass parameters, which are known only with large 
uncertainties. For this reason our estimations at this stage cannot give a preference to one 
or another of the above-mentioned hts to experimental data. 


E. Hadronic contributions to Aa 


Consider the hadronic correction to the electromagnetic hne structure constant a at the 
Z-boson scale. The evolution of the running electromagnetic coupling is described by 

a(0) 


a(s) = 


1 - Aaiept(s) - Aa[f2i(s) - Aah°i;(s) 


( 21 ) 


The leptonic part Aaiept(s) is known to the three loop level, Aaiept(M|) = 0.03149769 j^. 
It is conventional to separate the contribution Aa[|^^^(s) coming from the hrst hve quark 
flavors. The contribution of the f-quark is estimated as Aa|(°j(M|) = —0.000070(05) Q 
The quantity Aq;[|^j(s) at the Z-boson scale can be represented in the form of the dis¬ 
persion integral ©• The total function R{s) is 




( 22 ) 


where Q/ is the quark electric charge of flavour /. For the calculation of R{s) we use fl22jl 
with hve quark havors f = u,d, s, c, b. Varying the parameters as has been described above 
and using rric = 1.3-1.5 GeV, we get 


Mad(^l) = (278.2 ± 3.5) X 10 


-4 


(23) 


This value is to be compared with predictions extracted from a wide range of data describing 
hadrons l5lll: 


A«lad(^z) = (275.5 ± l.Oexpt ± 1.3rad) X 10 


-4 


(24) 


We see that our result (El is consistent with previous theoretical/experimental evaluations, 
with comparable uncertainties. 
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The relative error in (El is substantially less than the error that appears in the quantity 
and therefore one can obtain a more exact result. In comparison with the result, 
where the contribution of the c-quark was about 2 %, now it is about 30 %. The contribution 
of the 6-quark is about 5 % and the relative contribution of the t-quark is a fraction of a 
percent. 

IV. CONCLUSIONS 

A method of performing QCD calculations in the nonperturbative domain has been de¬ 
veloped. This method is based on the analytic approach to QCD, in which there are no 
unphysical singularities, and takes into account the summation of threshold singularities 
and the involvement of nonperturbative light quark masses. 

The following quantities have been analysed: the inclusive r-decay characteristic in the 
vector channel, the light-quark Adler function, D(Q^); the smeared i?A-function; the 
hadronic contribution to the anomalous magnetic moment of the muon, and the 

hadronic contribution to the fine structure constant, Aq;®^(M|). We have demonstrated 
that the proposed method allows us to describe these quantities rather well. 
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